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ABSTRACT 

A specialized data mining algorithm has been developed using wide- field photometry catalogues, en- 
abling systematic and efficient searches for resolved, extremely low surface brightness satellite galaxies 
in the halo of the Milky Way (MW). Tested and calibrated with the Sloan Digital Sky Survey Data 
Release 6 (SDSS-DR6) we recover all fifteen MW satellites recently detected in SDSS, six known 
MW/Local Group dSphs in the SDSS footprint, and 19 previously known globular and open clusters. 
In addition, 30 point source overdensities have been found that correspond to no cataloged objects. 
The detection efficiencies of the algorithm have been carefully quantified by simulating more than 
three million model satellites embedded in star fields typical of those observed in SDSS, covering a 
wide range of parameters including galaxy distance, scale-length, luminosity, and Galactic latitude. 
We present several parameterizations of these detection limits to facilitate comparison between the 
observed Milky Way satellite population and predictions. We find that all known satellites would be 
detected with > 90% efficiency over all latitudes spanned by DR6 and that the MW satellite census 
within DR6 is complete to a magnitude limit of My ~ —6.5 and a distance of 300 kpc. Assuming all 
existing MW satellites contain an appreciable old stellar population and have sizes and luminosities 
comparable to currently known companions, we predict lower and upper limit totals of 52 and 340 
Milky Way dwarf satellites within ~ 260 kpc if they are uniformly distributed across the sky. This 
result implies that many MW satellites still remain undetected. Identifying and studying these elusive 
satellites in future survey data will be fundamental to test the dark matter distribution on kpc scales. 
Subject headings: dark matter, dwarf galaxies. Local Group 



1. INTRODUCTION 

The dwarf galaxy population of the Milky Way (MW) 
provides invaluable insight into galaxy formation and 
evolution. Their resolved stars reveal their formation 
histories and enable precise measurements of their struc- 
tural parameters, ages and metallicities. These histories 
of individual, nearby systems provide a unique approach 
to studying the Universe across the cosmic ages. Dwarf 
galaxies are also the most numerous type of galaxy in 
the Universe and are thought to be the building blocks 
of larger galaxies. Owing to their low masses, their prop- 
erties may be strongly influenced by ionizing radiation in 
the early Universe and by the energy released by super- 
novae. The impacts of both of these are weak links in 
our understanding of structure formation. Finding and 
studying nearby dwarfs of the lowest masses and lumi- 
nosities is thus an essential component to understanding 
galaxy formation on all scales. 

The Milky Way dwarf galaxy population is also at 
present the most direct tracer of the abundance, mass 
spectrum, characteristic size, and spatial distribution of 
dark matter on sub-galactic scales. Standard ACDM 
simulations of MW-size dark matter haloes predict many 
more dar k matter sub-haloes than are observed as dwarf 
galaxies (jKlvpin et al.lll999HMo"ore et al.lll999 f). The re- 
cent "Via Lactea" simulation contains 2,000 dark matter 
sub- hales within 289 k pc of the simulated primary galaxy 
(jDiemand et al.l [20071 ) which have no observed optically 
luminous counterparts. This discrepancy gives rise to the 

^ Research School of Astronomy and Astrophysics, AustraUan 
National University, Cotter Road, Weston, ACT 2611; 

^ Clay Fellow, Harvard-Smithsonian Center for Astrophysics, 60 
Garden Street Cambridge, MA 02138 



questions how and in what mass regime do the baryons 
disappear from dark matter clumps? Studies of the spa- 
tial distributions of Milky Way and M31 dwarf galaxy 
companions have also highlighted possible discrepancies 
between ACDM theory and observations (jKroupa et al.l 
l200llKan7eraI1[2005tlMetz et al.l[200l . 

The most obvious reason for these apparent discrep- 
ancies in the number and spatial distributions of dwarf 
galaxies is substantial incompleteness as the Milky Way 
halo has not yet been uniformly searched for dwarf 
galaxy companions to low enough luminosities and star 
densities, in particular close to the Galactic plane 
where foreground con tamination is severe. For example, 
IWillman et al.l ()2004D compared the spatial distribution 
of MW satellites to that of M31's population, as well as 
that of a simulated dark matter halo and concluded that 
some dwarfs may have been missed at low Galactic lati- 
tudes, and that the total number of MW satellites with 
properties similar to the known objects could be as many 
as triple the known population. 

The viability of this solution, at least in part, has 
been underscored by the recent discoveries of four- 
teen new Galactic companions from the photometric 
data of the Sloan Digital Sky Survey (SDSS). These 
objects all appear to be dominated by old (> 10 
Gyr) stellar populat ions, with the exception of Leo T 
(|de Jong et al.l I2008D . Nine of these companions were 
immediately identified as dwarf spheroidal (dSph) galax- 
ies: Ursa Major, Canes Venatici, Bootes, Ursa Major 
II, Canes Venatici II, Hercul es, Leo IV, Coma Berenices, 
Leo T and Leo V (WiUmai reraI|l2005bl: [Zucker et alj 
'2006b'; 'Belokurov et al. 2006b'; Zucker et al' '2006d: 
i Belok urov et al. 2007; Irwin et al. 2007; Belokurov et al] 
l2008f ). Spectroscopic follow-up has confirmed that they 



all are highly dark mat t er dominated dwa rf galax- 
ies (iSimon fc Cxehal [20071 : iMartin et all [200l Will- 
man 1, Segue 1 and Bo5tes II (IWillman et al.1 l2005al : 
iBelokurovet al.ll2007t IWalsh et al.ll2007al) occupy a re- 
gion of size-luminosity space not previously known to be 
inhabited by old stellar populations, at the intersection of 
MW dSphs and glo bular clusters. Spectroscopic studies 
(jMartin et al.ll2007r ) showed that Willman 1 may be res- 
ident inside a dark matter subhalo with a mass-to-light 
ratio of ~ 470. If these ambiguous objects are gravi- 
tationally bound, then tidal a rguments also favor them 
being dark matter dominated (|Walsh et al.ll2007"bh . The 
remain ing two objects discov ered in SDSS, Koposov 1 
and 2, (jKoposov et al.|[2007b[ ) are extremely faint Galac- 
tic globular clusters. 

Numerous authors have shown that the predictions 
of ACDM simulations can be reconciled with the small 
number of observed MW dwarf galaxies if simple mod- 
els for baryonic physical processes are taken into ac- 
count when interp reting the resu l ts of numerical simu- 
lations (e.g. Bullock et all 120011 : iKravtsov erall |2004 
ISimon fc Gehall2007[ ). For example. IStrigari et 111(120071 ) 
show that the central masses (Mo.efepc) of the Milky Way 
dwarf galaxies are well-constrained by the data and that 
their mass function closely matches the Mo.ekpc rnass 
function of both the earliest forming sub-halos and the 
most massive accreted sub-halos in the Via Lactea sim- 
ulation. 

A well-defined, deep survey of Milky Way dwarf galax- 
ies over a large fraction of the sky is critical to assess any 
of the above scenarios. The dwarf galaxies detected (or 
not) by such a survey will provide one of the best ways 
to rigorously test the ACDM paradigm by comparing 
a variety of metrics (distribution, mass, scale sizes, and 
number) of the Milky Way dwarfs with the predictions of 
ACDM + galaxy formation models. Willman et al (2002) 
and iKoposov et al.l (J200 7a) have previously conducted 
automated searches for Milky Way dwarfs in the SDSS 
and their corresponding detection limits. The original 
Willman et al survey was only performed ov er a couple 
of thou sand square degrees of sky. The Ko posov et al.l 
(|2007ar ) survey was performed with a more sensitive al- 
gorithm (critical, because they found many new satellites 
to be on the edge of detectability), but few galaxies were 
used to accurately quantify their detection limits. 

In this paper, we present critical improvements to the 
present characterization of the detectability of Milky 
Way dwarf galaxies over the ~9,500 square degrees of 
the SDSS in Data Release 6 (DR6). We also present an 
improved detection algorithm over previous searches. We 
aim to construct the most complete, well-defined census 
of MW satellites by embarking on a Milky Way all sky 
satellite hunt. This search will ultimately combi ne SDSS 
Data Release 6 (jAdelman-McCarthv et"aIll2008D . The 2- 
Micron All Skv Survev (2MASS: ISbutskie et all [200^ 
and the upcoming Southern Sky Survey ([Keller et al.l 
l2007t l. The Southern Sky Survey will cover the entire 
20,000 deg^ below <5 < 0° using the new Australian Na- 
tional University (ANU) SkyMapper Telescope equipped 
with a 5.7 sq degree wide- field camera that is currently 
under construction with survey operation expected to 
commence early 2009. 

2. SDSS DATA 



The Sloan Digital Sky Survey (jYork et al.l l2000f) is 
an automated multi-color imaging and spectroscopic 
survey spanning 9,500 square degrees surrounding the 
North Galactic Pole. The u, g, r,i and z imaging data 
(jFukugita et al. "1996: G unn et 3111199 8) are photometri- 
cally an d astrometrically reduced through an automatic 
pipeli ne (JHogg et al.ll2001l : |Ivezic et al. 200l lSmith et all 
I2002t iTucker et alTl2006l : IPier'et al.i aOo£. We subse- 
quent ly correct for reddening with the ISchlegel et aP 
(1998) extinction values given in the SDSS catalog. 
All following work is performed on point sources from 
DR6, using the photometry flags from the examples 
of database queries appropriate for point sources avail- 
able on the SDSS Skyserver website'^. To ameliorate ef- 
fects of incompleteness in the point source catalog and 
star/galaxy separation, we only consider sources brighter 
than r = 22.0. The photometric data are provid ed by 
the SDSS DR6 (jAdelman-McCarthv et all [200l . We 
take our data from local copies of the SDSS datasets, 
maintained at the Harvard-Smithsonian Center for As- 
trophysics. 

3. SURVEY METHOD 

Low surface brightness Milky Way satellites are de- 
tectable only by their resolved stars. With the least lu- 
minous known Milky Way satellites, such as Bootes II, 
containing fewer th an ^ 20 stars bright er than the main 
sequence turn-off ([Walsh et al.l I2007H) . a deep, wide- 
field, uniform, multi-color photometric catalog is essen- 
tial for searching for these objects. They will typically 
reveal their presence as statistically significant spatial 
overdensities relative to the Galactic foreground. Their 
signal can be enhanced by selecting stellar sources that 
are consistent in color-magnitude space with, for exam- 
ple, an old population of stars at a fixed distance. In 
this paper, we restrict ourselves to the old stellar pop- 
ulations characteristic of Local Group dSphs, but the 
population-specific elements of the algorithm can be eas- 
ily modified for other systems. The stra tegy of our de- 
tection algorithm is built upon that of I Willman et al.l 
([2002) and Willmaa (2.003) which utilized the photo- 
metric catalogs fr om SDSS and led t o the discoveries 
of Ursa Major (Willm an et all [20Ci5bl) and Willman 1 
(IWillman et al.l 2005al. I t is also similar in sp irit to 
iBelokurov et all (|2007i ) and lKoposov et all ([2007a[ ). Sev- 
eral systematic searches for MW dwarfs h ave also been 
done with non-SPSS da ta ([IrwinI [19941 : iKlevna et all 
[19971 : IWhiTing et al.l[200l . 

In summary, our algorithm applies color and magni- 
tude cuts to stars in the DR6 catalog, stores their distri- 
bution in a spatial array with 0.02° x 0.02° pixels, spa- 
tially smooths the array with a Plummer surface density 
profile, and sets comprehensive thresholds for detection. 
Each of these steps is described in detail in the following 
sections. 

3.1. Data Management 

In order to efficiently manage thousands of square de- 
grees of survey data in a catalog containing tens of mil- 
lions of stars, we ffist divide the dataset (in the case 
discussed in this paper, SDSS DR6) into stripes, each 
spanning 3° in declination (to avoid projection effects) 

^ http://cas.sdss.org/dr6/ 



with 2° of overlap in Dec between adjacent stripes. This 
overlap creates a substantial redundancy to ensure that 
real objects are situated in the central ~ 2° of Declina- 
tion in at least one stripe, away from possible edge-effects 
introduced at the stripe boundaries during the processing 
described in §3.3. We then take the longest continuous 
regions of the DR6 footprint in RA. 

3.2. Selection Criteria 

The mainly old, metal-poor stars of a nearby dwarf 
galaxy will occupy a well defined locus in the color- 
magnitude diagram (CMD), in contrast to Milky Way 
stars which span a wide range in distance, age, and 
metallicity. Therefore selecting stars that are consistent 
in color-magnitude space with a population of old stars at 
a particular distance will significantly enhance the clus- 
tering contrast of a dwarf galaxy's stars over the fore- 
ground noise from Milky Way stars. 

We use th e oretic al isochrones in SDSS filters from 
iGirardi et al.l ()2004f ) to define the regions of (r, g ~ r) 
space likely to be po pulated by old, metal-poor stars. 
ISimon fc Gehal (|2007l ) obtained spectra of stars in eight 
of the newly discovered dwarfs: CVn, CVn II, Com, Her, 
Leo IV, Leo T, UMa and UMa II and found mean abun- 
dances in the range —2.29 <[Fe/H]< —1.97. Based on 
this result, we consider isochrones for populations with 
a bundances of [Fe/H ] = — 1.5 and —2.27 (the lower limit 
in IGirardi et al . 2004) and with ages 8 and 14 Gyr. Four 
isochrones in these ranges can be used to bound the re- 
gion of CMD space we are interested in, namely the four 
combinations of [Fe/H] = —1.5 and —2.27 and ages 8 and 
14 Gyr. Figure [l] shows these four isochrones projected 
to a distance of 20 kpc. 

We define the selection criteria by the CMD envelope 
inclusive of these isochrones +/- the 1 sigma (g — r) color 
measurement error as a function of r magnitude. Shift- 
ing these isochrones over distances between m — M = 
16.5 - 24.0 in 0.5 mag steps defines 16 different selection 
criteria appropriate for old stellar populations between 
d ^ 20 kpc and > 630 kpc. We truncate our color- 
magnitude selection template at a faint magnitude limit 
of r = 22.0, beyond which photometric uncertainties in 
the colors and star/galaxy separation limit the ability to 
detect these populations. We also truncate the selection 
template at g — r = 1.0 as including redder objects adds 
more noise from Milky Way dwarf stars than signal from 
more distant red giant branch (RGB) stars. Finally we 
do not include stars with Sg or Sr > 0.3 mag in our analy- 
sis. To efficiently select stars within this CMD envelope, 
we treat the CMD as an image of 0.025 x 0.125 (color 
X mag) pixels and determine which stars fall into pixels 
classified as "good" according to the selection criteria. 
Figure [T] shows an example of the selection criteria, in 
this case for m — M — 16.5 (~ 20 kpc). The shaded 
region highlights pixels that would be classed as "good" 
for a system at ~ 20 kpc. 

3.3. Spatial Smoothing 

After the photometric cuts are applied, we bin the 
spatial positions of the selected stars into an array, E, 
with 0.02° X 0.02° pixel size. We then convolve this 
2D array with a spatial kernel corresponding to the ex- 
pected surface density profile of a dSph. We refer to 
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Fig. 1. — (r, g — r) color-magnitude diagram showing the two 
reddest and two bluest theoretical isochrones for old stellar pop- 
ulations ([Fe/H]= —2.27, —1.5 and age 8, 14 Gyr) at a distance 
modulus oi m — M = 16.5 {^ 20 kpc), generated from Girardi et 
al. (2004). The shaded region shows pixels that pass the selection 
criteria. 



TABLE 1 
Angular sizes of the 
satellites detected in 

SDSS. 



Object 


rh 




(arcmin) 


Bootes 


12.6 


Bootes II 


4.2 


Canes Venatici 


8.9 


Canes Venatici II 


1.6 


Coma Berenices 


6.0 


Hercules 


8.6 


Leo IV 


2.5 


Leo V 


0.8 


LeoT 


1.4 


Segue 1 


4.4 


Ursa Major 


11.3 


Ursa Major II 


16.0 


Willman 1 


2.3 



this smoothed spatial array as A. For our spatial ker- 
nel we use a Plummer profile with a 4.5' scale-length. 
This value provides an effective compromise between the 
angular scale-lengths of compact and/or distant objects 
with those of closer/more extended objects. For refer- 
ence the angular sizes of the new satellite s are listed in 
Table [H We use the rh values derived by iMartin et al.l 
(|2008[ ) except for Leo V (jBelokurov et al.ll2008f ). 

The normalized signal in each pixel of A, denoted by S, 
gives the number of standard deviations above the local 
mean for each element: 



S^ 



A- A 



The arrays of running means, A, and running standard 



deviations, A^r are both calculated over a 0.9° x 0.9° win- 
dow around each pixel of A. In particular, Ag- is given 
by: 



A„ 



ln{A-Ay*B-{{A-A)*B) 



n(n — 1) 



i? is a box-filter with n elements and is the same size 
as the running average window. The resulting array A^^ 
gives the standard deviation value for each pixel of A 
as measured over the 0.9° x 0.9° span of the filter. In 
the next section, we will define the detection threshold 
of this survey in terms of 5, as well as in terms of the 
local stellar density E. 

3.4. Detection Threshold(s) 

In a large survey such as ours, it is critical to set detec- 
tion thresholds strict enough to eliminate false detections 
but loose enough to retain known objects and promising 
candidates. To characterize the frequency and magni- 
tude of truly random fluctuations in stellar density an- 
alyzed with our algorithm, we measure the maximum 
value of S for 199,000 5.5° x 3° simulated fields of ran- 
domly distributed stars that have been smoothed as de- 
scribed in the previous section. The only difference is 
that there is no gradient in stellar density across each 
field. In the interest of computational efficiency we do 
not use a running window for the mean and a of each 
simulated field. The field size is chosen such that 1,000 
fields roughly totals an area equal to the DR6 footprint 
(neglecting regions lost during convolution). We select 
199 stellar densities ^n.^, to simulate, linearly spaced be- 
tween 10 and 4,000 stars per square degree. This range 
of stellar densities is to model the density range we find 
after applying the color-magnitude selection criteria de- 
scribed in ii3.2l across the SDSS. In i j5.1l we study the 
variation of detection limits with Galactic latitude (fore- 
ground stellar density); the typical number densities we 
will consider there are higher than 10 - 4,000 stars per 
square degree because we wish to parametrize the detec- 
tion limits in terms of the density of all stars bluer than 
g — r = 1.0 and brighter than r = 22. 

Figure [2] shows a two dimensional cumulative his- 
togram of the 199,000 max(S') values as a function of 
stellar density. In low density fields, the distribution of 
pixel values becomes non-Gaussian so a simple, global 
threshold value is insufficient. The solid grey line shows 
the contour containing 99% of the 199,000 max(5') values 
at each density. If we simply used a value like this as our 
threshold, we would be biasing ourselves against detect- 
ing extended objects: Large angular scale-length systems 
may not have a peak pixel value above this value, for ex- 
ample because stars in the object itself increases the local 
running mean and sigma. However, such an overdensity 
may have a characteristic area larger than any random 
fiuctuation. 

We thus define a detection threshold based on both the 
peak density and the characteristic area of an overden- 
sity. To define such an area we scale down the 99% con- 
tour from Figure[2]and define a threshold density Sth{n*) 
as a function of stellar number density. Then, using the 
199,000 random fields we examine the relationship be- 
tween the peak density max(S') divided by Sth{n*) and 
"detection" area, i.e. the area of contiguous pixels of S 



that have values above Sthin*) (the white line). Figure 
[3] shows this area versus max(S')/5't/i(ri,) for the ran- 
dom fields. If we assume a purely random foreground, 
we would expect one false detection in the DR6 foot- 
print above an area of ~ 55 square arcminutes or above 
a peak density of ~ 1.6 x Sth{n*)- These numbers are 
set by the factor by which we scale down the threshold 
function and are themselves arbitrary. 




Fig. 2. — A two dimensional cumulative histogram showing 
the distribution of max(S') values for smoothed fields as a func- 
tion of different stellar densities n,. The grey line bounds 99% of 
the niax(5) values and the white line shows our threshold density, 
Sth = f{n*). 
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Fig. 3. — "Detection" area versus niax{S')/5t;i{n,) for the 
199,000 random fields. The black contour shows the level at which 
purely random clustering would produce one false detection over 
the approximate area of DR6. 

Based on the results of these simulations we set the 
area threshold to a more conservative 60.0 sq arcmin 
and the density threshold to a more conservative 1.75 x 
Sthin*) to eliminate false positive fiuctuations while pre- 
serving all of the known obje cts within DR6, including 
Bootes II (IWalsh et al.ll2007af) and the Koposov 1 and 2 
globular clusters (jKoposov et al.ll2007bD . Thus, a detec- 
tion is defined as a region where: 



• the area of a group of pixels of S contiguously 
above Sthin*) (white hue, Fig [2|) is greater than 
60.0 square arcminutes 

or 

• any single pixel value is greater than 1.75 x Sthin^,)- 

We implement these adaptive density thresholds as a 
function of local stellar density n* , so that the algorithm 
may be run over large fields with varying density and 
allow direct comparison between fields of greatly different 
densities. The stellar density n* is calculated for each 
pixel of the smoothed, normalized, spatial array S, as 
the 0.9° X 0.9° running average of the original spatial 
density array E. 

To summarize our algorithm: 

• Apply CMD cuts, bin spatial positions of remaining 
stars into E 

• Smooth E with Plummer profile to get A 

• Calculate the 0.9° x 0.9° running mean A and run- 
ning standard deviation Aa- 

• Define S" as S" = (^ - A)/A^ 

• Calculate array of threshold values Sth as function 
of stellar density n* (from 0.9° x 0.9° running mean 

of £;) 

• Detections are where contiguous regions of pixels 
with S > Sth is greater than 60.0 sq arcmin or any 
single pixel is greater than 1.75 x Sth- 

3.5. Identifying and Evaluating Detections 

For each of our DR6 data strips defined in §3.1, the 
steps outlined in the previous sections are repeated in 
0.5 magnitude distance modulus intervals, and these 16 
frames are layered to form a 3-dimensional array. This 
3D approach eliminates complications with multiple de- 
tections of a single object using selection criteria for dif- 
ferent distance moduli, and selects out the strongest de- 
tection. The coordinates of stars within each detection 
and the CMD within the detection's area are plotted for 
later visual inspection. Galaxy clusters and point sources 
around partially resolved background galaxies (such as 
their associated globular clusters) will contaminate the 
detections, but these can be identifiable based on their 
CMDs (see FigureOin g]), leaving a list of potential new 
Milky Way satellite galaxies and globular clusters. At 
this point follow up observations are typically necessary 
to confirm the existence and nature of these candidates. 

4. APPLICATION TO SDSS DATA RELEASE 6 

We apply our search algorithm (as described in §3) to 
21,439,777 sources with r < 22.0 and g - r < 1.0 in the 
9,500 deg^ of imaging data in Data Release 6 of the SDSS. 
The DR6 footprint is shown in Figure El along with pre- 
viously known dSphs (open blue circles) and satellites 
discovered in SDSS (closed red circles). 

The significance of our detections of known objects in 
terms of their peak density and area are shown in Figure 
[SI In the total area of DR6 analyzed, we find 100 unique 
detections above the thresholds, defined by the dotted 
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Fig. 4. — Aitoff projection of the DR6 footprint in Galactic 
coordinates, centered on the Galactic center. Previously known 
dwarfs are marked with open blue circles, satellites discovered in 
SDSS are marked with filled red circles. 



lines of Figure \5[ The positions of each of these detec- 
tions are cross-referenced against the SIMBAD database 
* as well as visually inspected via the SDSS Finding 
Chart Tool^. Of our 100 detections, 19 are MW/Local 
Group dwarfs (counting Bootes II, Willman 1 and Segue 
1), 17 are Galactic globular clusters (including Koposov 
1 and 2), 2 are known open clusters, 28 are clustering of 
point sources associated with background galaxies such 
as unresolved distant globular clusters, and four are Abell 
galaxy clusters. The remaining 30 do not correspond to 
any catalogued objects, but color-magnitude diagrams of 
only a handful of these are consistent enough with a faint 
MW satellite to warrant follow-up. The remainder may 
be galaxy clusters whose detected center differs from its 
cataloged centre by more than ~ 0.25°, or perhaps tidal 
debris. If the MW stellar halo is the result of accretion 
of dSph then evidence of this accretion is expected. It 
should be noted that objects with relatively large angu- 
lar size, such as Draco and Sextans, substantially increase 
the average stellar density of the area they occupy which 
increases the threshold density, meaning they are not as 
high above the density threshold as one might expect. 
Due to the area threshold however, they are still very 
prominent detections. 

We recover all of the newly discovered objects that are 
within DR6 and the "classically" known Draco, Leo, Leo 
II, Leo A, Sextans, and Pegasus DIG dwarfs. Our detec- 
tions of the new dwarfs are presented in Figures El H and 
[8l These figures arc identical to those output by the auto- 
mated algorithm for each detection, aside f rom the addi- 
tion o f figure titles (My and distances from Martin et ahl 
|2008| and references therein). The left panel shows the 
spatial positions of stars passing the photometric selec- 
tion criteria at the distance modulus the object was most 
strongly detected at. The middle-left panel shows the 
contour plot corresponding to S, where the contour lev- 
els are iS)/Sth{n*) = 1.0, 1.2, 1.4, 1.6, 1.8, and 2.0. 

"^ http://simbad.u-strasbg.fr/simbad/ 

^ http://cas.sdss.org/astrodr6/en/tools/chart/chart.asp 
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TABLE 2 
Positions of Strongest MW Satellite Candidates 



Designation 


a 


<5 


(a, 6) 


D (kpc) 


"CVn W" 
"Her X" 
"UMa Y" 
"Vir Z" 


13:16:04.8 
16:27:45.6 
12:11:57.6 
12:20:19.2 


+33:15:00 

+29:27:00 

+53:35:24 

-1:21:00 


(199.02,33.25) 
(246.94, 29.45) 
(182.99,53.59) 
(185.08,-1.35) 


~ 160 

~ 100 

~ 100 

~40 



Fig. 5. — Same as Figure[3]but showing all detections in DR6. Dotted black lines show the adopted thresholds. Galactic/Local Group 
dSphs and Koposov 1 and 2 are shown as black filled circles. The brightest objects such as Draco and Sextans increase the average stellar 
density of the area they occupy, which increases Sth(j^*) in their vicinity. This means that max(S')/5'th('^*) is not as high as one might 
expect. 

The middle-right panel is the CMD of the detection area 
and the right panel is the field subtracted Hess diagram. 
The isochron e is th at of a 13 Gyr, [Fe/H]= —2.3 from 
iGirardi et al.1 ()2004D at the distance specified. Besides 
demonstrating the effectiveness of our algorithm, these 
detections provide a benchmark with which to compare 
candidates and determine which are consistent with be- 
ing a new dwarf satellite. 

To further illustrate the product of our algorithm we 
also show examples of undesired detections in Figure [51 
the galaxy cluster Abell 1413 (top) and Virgo cluster 
galaxy NGC 4486 (bottom). These represent typical de- 
tections of background galaxies and galaxy clusters. 

4.1. Candidate Milky Way Satellites 

Figure [10] shows four unidentified overdensities that 
have CMDs qualitatively similar to that of a dSph. 
All show statistically significant spatial clustering and 
do not coincide with a visible overdensity of back- 
ground galaxies. While several unknown detections have 
CMDs broadly consistent with old stellar populations, 
we present here four detections that are as strong or 
stronger than the detections of UMa and PegDIG. We 
present their positions and distances estimated from the 
CMD in Table [H 

5. EXPLORING DETECTION EFFICIENCY WITH 
SYNTHETIC SATELLITES 

The most advantageous aspect of a large, uniform 
search for MW dwarfs is the ability to rigorously cal- 
culate its detection limits in order to compare observa- 
tions with predictions. To calculate the detection com- 
pleteness of our search, artificially generated galaxies are 



embedded in simulated stellar foreground fields and put 
through the detection algorithm to investigate the sensi- 
tivity as a function of galaxy distance, luminosity, scale- 
length, and Galactic latitude. In this section, we describe 
in detail the method used to synthesize artificial SDSS 
fields and dSph satellites. 

5.1. Sowing the Simulated Fields 

When simulating fields in which to embed artificial 
galaxies, our goal is to create a large number of fields 
with the same point source color, magnitude, and den- 
sity distributions as observed in the SDSS DR6 footprint. 
The detectability of a dSph may change depending on its 
position in the sky. For example, those at low Galactic 
latitudes will be harder to detect than those at high lat- 
itudes, owing to the greater number of foreground stars. 
The relative proportions of the thin disk, thick disk, and 
stellar halo will also vary with latitude and longitude, 
changing what fraction of foreground stars will be in- 
cluded in the color-magnitude selection described in §3. 

To conduct a controlled experiment to see how Galac- 
tic foreground affects detection efficiency over the DR6 
footprint, we first select three fiducial latitudes to simu- 
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Fig. 6. — Our detections of recently discovered MW satellites. Left: Spatial plot of sources passing selection cut. Middle Left: Contour 
of smoothed spatial plot. Contours show 1.0, 1.2, 1.4, 1.6, 1.8, and 2.0 times the density threshold. Middle Right: CMD of region enclosed 
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Fig. 7. — Detections of recently discovered MW satellites cont. 
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latitudes above and below which 10% of the survey lies. this fraction by cos 6 gives the relative area on the ce- 

Figure [Tl] shows the fraction of sky observed by SDSS lestial sphere that each observed latitude occupies (solid 
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Fig. 9. — Detections of Abell 1413 (top) and NGC 4486 (bottom) as examples of galaxy cluster and background galaxy detections. 



grey line), showing that the majority of the DR6 foot- 
print by area is located between h « 45° and b « 65°. 
The cumulative total (solid black line) allows us to choose 
the latitudes corresponding to 10, 50 and 90% levels of 
DR6, namely 31°, 53° and 74°. These are the three val- 
ues of latitude that we implement in our simulations. 

Now that we have chosen what latitudes to simulate, 
we need to relate these to the stellar foreground den- 
sity. Figure [12] presents a 2D histogram of latitude and 
foreground density, considering only stars brighter than 
r = 22.0 and bluer than g — r = 1.0. This figure shows 
a span in foreground levels at each latitude. The solid 
black line traces the median and our chosen latitudes are 
marked along the x-axis. For each of our latitudes, we 
take a slice through the 2D histogram and use this dis- 
tribution of densities to randomly assign a density for 
each of our simulated fields. Each artificial star in our 
simulated fields is assigned photometric parameters from 
a star in DR6, chosen at random from all stars within 
±0.5° of the latitude in question. These stars are then 
randomly distributed in a 3° x 3° field. 

5.2. Forging Virtual Dwarfs 

To simulate a dSph galaxy CMD, we enlist HST obser- 
vations of thre e MW sat ellites^ : Carina, Draco, and Ursa 
Minor (Ho ltzman et al.li2006) . Figure [T51 shows the com- 
bined CMD of these objects. We take this My and V — I 
CMD and tr anslate it in c o lor an d magnitude loosely 
to match the iGirardi et al.l (|2004f ) isochrones g and r. 
The reasonable agreement betwee n the HST data af - 
ter approximate transformation and lGirardi et al] (|2004l ) 
isochrones in g and r is sufficient to use the HST data 
for our simulated objects. 

We use these data to create a composite old stellar pop- 
ulation catalog of stars brighter than Mr = 6 by combin- 
ing sources from the three HST dwarfs. Carina, Draco, 
and Ursa Minor each contribute 5,548, 4,487, and 3,296 

® http://astronomy.nmsu.edu/holtz/archival/litml/lg.html 



stars respectively. Each time we simulate a dwarf galaxy 
of n stars, we select those n stars at random from this 
composite catalog. The luminosity is calculated from the 
integrated flux from all stars, and a correction added to 
account for stars below an absolute magnitude of Mr = 6. 
The cumulative luminosity functions in the right panel 
of Figure \W\ show that typically ~ 10% of the total flux 
originates from stars below this cut-off. We then adjust 
the photometry of the stars to the correct distance mod- 
ulus and add photometric scatter to reflect increasing 
measurement uncertainty with fainter magnitudes. We 
do this by finding the best fit for the la magnitude un- 
certainty as a function of magnitude for each band in the 
SDSS data, and adding a normally distributed random 
realization S of this value cr(r7i) for the adjusted magni- 
tude, rristar = m + Sa{m). We then assign random posi- 
tions based on a Plummer surface brightness profile with 
a specified physical scale length at a given heliocentric 
distance. 

Figure [14] shows examples of three simulated dSphs. 
The middle panel shows a system not unlike Bootes H, 
highlighting the paucity of stars in the objects we are 
searching for. It is important to note that at these low 
luminosities, the total luminosity of galaxies with the 
same number of stars can vary dramatically, with this 
variation increasing for galaxies with fewer stars. A sin- 
gle RGB star can have a magnitude of Mr = —3 which 
is well in the regime of the total magnitudes of recently 
discovered satellites. Each generated galaxy is embedded 
in a simulated field, and then processed as described in 
§3.2, §3.3 and §3.4. 

5.3. Charting Detection Efficiency 

To test the efficiency of our search algorithm as a 
function of galaxy luminosity, scale-length, distance, 
and Galactic latitude, we generate a total of 3,825,000 
galaxies spanning ranges in Galactic latitude, luminosity, 
physical size, and distance. We simulate systems at lati- 
tudes of 31°, 53° and 73° and with 2^ x 100 stars brighter 
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Fig. 10. — The four unknown detections with detection strength equal to or greater than known satellites (except Leo V). From top 
to bottom: "Canes Venatici W", "Hercules X", "Ursa Major Y" and "Virgo Z". Isochroncs show the distance interval at which these 
overdensities produced the strongest detections. 

than My = 6, where x is an integer between and 11 
(giving a range of 100 to 204,800 stars). These stellar to- 
tals correspond to mean total magnitudes of My = —1.5, 
-2.3, -3.1, -3.9, -4.7, -5.5, -6.2, -7.0, -7.7, -8.5, 
—9.2 and —10.0. For each of the 36 combinations of 
latitude and magnitude, we simulate a large number of 
galaxies with distances and physical scale-lengths ran- 
domly generated with the limits 1.3 < \ogd/kpc < 3.0 
and 0.9 < log r/j/pc < 3.0. For the brightest and faintest 
systems we tailor these ranges to avoid redundant it- 
erations; there is little to be gained by simulating a 
My = -1.5 system at 200 kpc or a My = -10.0 sys- 
tem at 20 kpc. Hence, the total number of simulations 
for each magnitude/latitude combination varies, but is 
chosen such that there are typically 500 simulations in 
each 0.1 log(d) x 0.1 log(r;i) bin of Figure [TSl 



6. DISSECTING EFFICIENCY TRENDS 



Figure [12] shows the detection efficiency of simulated 
dwarf galaxies as a function of luminosity, scale-length, 
and distance at the median SDSS Galactic latitude (53°). 
Overplotted are each of the Milky Way satellites detected 
in SDSS, not including Koposov 1 and 2. Each panel 
contains a grey-scale map of the detection efficiency for 
simulated galaxies of the mean absolute magnitude spec- 
ified in the panel. Because the total magnitude varies 
for systems with a constant number of stars, we quote 
both the mean magnitude and the standard deviation 
of magnitudes for each panel. White shows regions of 
100% efficiency, while black shows 0%. The four con- 
tours, moving outwards from 100% efficiency, show the 
90%, 84.13%, 50%, and 15.86% levels. The 84.13% and 
15.86% levels were chosen to illustrate the ilcr in de- 
tectability as a function of distance and size. 

The greater the number of stars in a simulated galaxy. 
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Fig. 11. — The fraction of sky as observed by SDSS DR6 at each 
latitude. The dotted line shows what fraction of the small circle 
on the celestial sphere traced by each latitude has been surveyed, 
and the grey line is this fraction weighted by the cosine of latitude, 
to give a relative sky area observed at each latitude. The largest 
area of DR6 observations occur at 6 ~ 60° . The solid black line is 
the cumulative total of the grey line. 
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Fig. 12. — Grey-scale plot of the number of sources per square 
degree satisfying r < 22.0 and g — r < 1.0 versus absolute Galactic 
latitude. The black line traces the mean density at each latitude. 
The three latitudes we choose to simulate are marked along the 
X-axis. 

the less its absolute magnitude will vary between realiza- 
tions, so the standard deviation in the integrated mag- 
nitudes of simulated galaxies contributing to each panel 
decreases with increasing magnitude. The My = —7.7 
panel which shows a small increase in standard deviation 
marks a change in the way the galaxies are simulated; we 
are now simulating systems with more stars than are in 
our HST catalog so the simulated stars no longer have 
unique photometry drawn from this catalog. This am- 
plifies the small number effect of single stars on the total 
magnitude. 

Figure [15] illustrates the necessity for a large number 
of simulations as there are subtle features that would 
otherwise be unresolved. For example, this figure shows 
that the detectability of dwarfs is not a step function in 
distance, but rather slowly falls off at a rate that differs 
for systems with different t otal lum inosities. This is in 
contrast to iKoposov et al.l (|2007af ). who found a steep 
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Fig. 13.— HST data of thr ee Milky Way satellite dSphs (Ca- 
rina, Drac o, and Ursa Minor ; IHoltzman et al] 120061 ) with SDSS 
isochrones jGirardi et al.||2004j ) overlaid. The right panel shows the 
cumulative luminosity functions for the corresponding isochrones, 
using the four combinations of [Fe/H]= —2.27,-1.5 and age= 
8, 14Gyr. Data are corrected for distance and presented in absolute 
magnitude. 

boundary between 100% and 0% efficiency. The gradual 
fall-off in dwarf detectability with distance will be dis- 
cussed in more detail in ij6.1l b u t we p ostpone detailed 
comparison with 'Kopo sov et al.1 (|2007a[ ) until ^3 

The critical factor affecting the detectability of an ob- 
ject with our algorithm is the number of stars brighter 
than r = 22.0 that fall under the rh = 4.5' Plummer 
smoothing kernel. In the following sections, we use this 
to gain physical understanding of the features of Figure 
1151 and to derive an analytic expression to describe de- 
tection efficiency as a function of galaxy magnitude, size, 
distance, and Galactic latitude. This analytic expression, 
as well as a routine to interpolate detectability directly 
from the simulations, will be made publicly available and 
can be used, for example, to correct the MW satellite lu- 
minosit y function as previously done by iKoposov et al.l 
()2007a| ) or make an estimate on the corrected radial dis- 
tribution of MW satellites. Such endeavours are beyond 
the scope of this paper, but we use the function to es- 
timate the total number of MW satellites that remain 
undetected, presented in fj8l 

6.1. Efficiency versus Distance 

Figure [15] shows that the detectability of resolved 
dwarfs around the Milky Way is not a step function 
in distance. As distance to a dwarf galaxy increases, 
the number of stars brighter than r = 22.0 {Nr<22) de- 
creases. In an idealized scenario, the detectability of that 
dwarf would drop from 100% to 0% at a distance beyond 
which the number of resolved stars required to produce a 
detection is larger than Nr<22- As we have discussed in 
previous sections, random variations in the stellar lumi- 
nosity function can be substantial in the faintest systems, 
hence Nr<22 will be affected by stochastic fluctuations. 
Moreover, the wide range in foreground densities at a 
given Galactic latitude (see Figure [T2|) impacts the de- 
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Fig. 14. — Simulated dSph systems. 
Top: d = 25 kpc, r^ = 250 pc, My = -4.7. 
Middle: d = 45 kpc, r^ = 40 pc, My = -2.3. 
Bottom: d = 100 kpc, r^ = 250 pc, Mv = -5.6. 

tectability of two identical dwarfs. Therefore the transi- 
tion from 1.0 to 0.0 detection efficiency is not expected 
to be a step f unction, but ra t her des cribed by a Gaus- 
sian integral. i Koposov et al.l (j2007a[ ) also modeled the 
detectability transition with a Gaussian integral, despite 
finding a steep decline. So detection efficiency as a func- 
tion of log(distance), which for brevity we denote ^d can 
be described as 



DE{Ld) = erf 



■4^ 



<^Ld 



where id is the logarithm of distance, and id and a^d 
are the mean and standard deviation of log(distance). 
The mean corresponds to the distance at which a sys- 
tem would be detected with 50% efficiency. The error 
function erf is defined as: 



erf{x) 



Vi^) 



dt. 



Examination of Figure [T51 shows that a^d (how quickly 
efficiency transitions from unity to zero) changes depend- 
ing on luminosity. We would naively expect a^a to con- 



tinue increasing with decreasing brightness as small num- 
ber statistics becomes more dominant. Instead it shows 
a maximum at My ~ —3.9 before decreasing. This is a 
result of the stochastic fluctuations and the derivative of 
the luminosity function; since Nr<c22 varies for systems 
with the same total number of stars, the individual dis- 
tance that each of these systems could be detected at also 
varies. As the number of stars above the brightness limit 
is dependent on the LF, the slope of the LF determines 
how Nr<:22 changes with distance. Hence a^d is smaller 
for the faintest objects when the MSTO is required for a 
detection because the LF is at its steepest at the turnoff. 

6.2. Efficiency versus Scale-length 

The fraction of a dwarf's stars within our 4.5' spa- 
tial smoothing kernel decreases with increasing physi- 
cal scale-length and/or decreasing distance. A system of 
some luminosity and distance that is detectable when its 
angular size is < 4.5' may thus be undetectable if those 
same stars are spread over a larger angular scale. As the 
concentration of stars increases we would expect detec- 
tion efficiency to also increase. However, when the angu- 
lar size of a dwarf is comparable to the smoothing kernel 
size, the detectability does not appreciably improve with 
further decrease in size since the number of stars within 
the Kernel is not significantly changing. Hence objects 
of this angular size or smaller will be detected with the 
same efficiency. 

Once the angular size becomes larger than the kernel 
size, the number of stars within the kernel declines. The 
relationship between size and detection distance is de- 
pendent on the stellar luminosity function of the system. 
Take for example an object with an angular size larger 
than the smoothing kernel, detected with 50% efficiency 
at some distance. To keep the object at 50% efficiency 
as we continue to increase the physical size, the drop in 
efficiency can be counteracted by decreasing the object's 
distance. As this object is moved closer the number of 
stars above r = 22.0 increases at a rate corresponding to 
the LF. At the distance when the main sequence turn- 
off (MSTO) becomes brighter than r = 22.0 (~ 65 kpc), 
the rapid increase in the number of stars corresponds to 
a sudden improvement in detection efficiency, evident in 
Figure [15] at log(d) « 1.8 in the My = -5.5, -4.7 and 
—3.9 panels. As with distance, efficiency versus scale- 
length can be modeled by a Gaussian integral, but with 
the mean and standard deviations as functions of dis- 
tance; so 



DE{Lrh) = erf 



(^LrALd) 



6.3. 



Analytically Expressing Detection Efficiency 

Combining the previous two results, we can analyti- 
cally describe detectability with a Gaussian integral over 
logd multiplied by another Gaussian integral over logr^. 
To introduce magnitude My and latitude h, we set the 
means (id, if/j) and standard deviations (cr^d, (^^ru) in 
the integrals to be functions of My and h. Therefore the 
detection efficiency DE can be expressed: 



DE = erf 



J{My,b)^Ld 

a^d{Mv,b) 



X erf 



'^rhJLdjMy ,b)) ~L Th 
a,r,ALd{My,b)) 
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Fig. 15. — Detection Efficiency for specific galaxy parameters. Each panel shows the detection efficiency as functions of distance and 
scale-length for a particular number of galaxy stars. The average total absolute magnitude for each of these sets of galaxies is shown 
along with the standard deviation in magnitudes. Contours show the 90%, 84.13%, 50%, and 15.86% levels. Sizes and distances of known 
MW dSphs are shown in the best matching magnitude panel. Leo and Le o A are shown in gr ey as they are sig nificantly brighter tha n 
My = —10.0. Values for newly discovered objects plus Draco are taken from lMartin et aTl II2008I ) except for Leo V IIBelokurov et ani2008l ). 
Values for Leo, Leo II, Leo A, and Sextans taken from iMateol llT998i ). 
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TABLE 3 

Comparison of Interpolated and 

Analytical Detection 

Efficiencies 



Object 


Intcrp. 


Analyt. 


Diff. 


Boo 


99.84 


99.98 


-0.14 


Boo II 


90.39 


95.96 


-5.57 


CVn 


100.0 


100.0 


0.0 


CVnII 


98.26 


99.07 


-0.81 


Com 


99.71 


100.0 


-0.29 


Dra 


100.0 


100.0 


0.0 


Her 


99.82 


99.83 


-0.01 


Leo 


100.0 


100.0 


0.0 


Leo II 


100.0 


100.0 


0.0 


Leo A 


100.0 


100.0 


0.0 


Leo IV 


98.65 


99.77 


-1.12 


Leo V 


83.56 


91.27 


-7.71 


Leo T 


93.41 


99.38 


-5.97 


Segue 1 


100.0 


100.0 


0.0 


Sex 


100.0 


100.0 


0.0 


UMa 


99.97 


99.86 


0.11 


UMall 


100.0 


99.96 


0.04 


Will 1 


98.56 


99.30 


-0.74 



The means and standard deviations can be found by 
fitting Gaussian integrals along the distance and scale- 
length axes of the panels in Figure [TSl The main source 
of uncertainty in this expression is the galaxy luminosity, 
which while correlated with the number of stars, can vary 
by over a magnitude for systems of equal detectability. 
The function does however give a good statistical approx- 
imation from which to estimate the properties of the true 
Milky Way satellite population. In FigurefTHlwe compare 
the My = —3.9 panel of Figure [T51 with the analytical 
function. There is good agreement between the empirical 
and analytical detection efficiencies with a la deviation 
of only ~ 8.7% across the entire range of parameters. 
For reference the size of the 4.5' smoothing kernel and 
the distance at which the MSTO becomes resolved are 
shown in the center panel in red and blue respectively. 

The analytical efficiency is compared in Table [3] to the 
interpolated efficiency from the grids in Figure fTSl for real 
MW dwarfs. For all objects besides Boo II, Leo V, and 
Leo T, the difference is within ~ 1%. At first inspection 
it may seem odd that all objects have a very high, al- 
most 100%, efficiency, but given that most of the param- 
eter space probed by our simulations yields either zero 
or unity efficiency, it is not unexpected that the handful 
of objects detected in this vast volume are detected with 
high efficiency. 

6.4. Efficiency versus Latitude 

Unlike the other parameters which vary a dSph's sig- 
nal strength. Galactic latitude affects detection efficiency 
by changing the foreground density, and therefore noise 
above which we must detect a signal. If latitude plays a 
significant role in the detectability of dwarfs, then it must 
be taken into account when making any corrections to the 
MW satelhte census. Figure [17] shows the My « -3.9 
panel of Figure [12] {b = 53°) with the addition of the 
50% detection efficiency (dashed) and 90% detection ef- 
ficiency (dotted) contours of the b — 31° (orange) and 
74° (blue) simulations ovcrplotted. As expected, an ob- 
ject of given size and luminosity will not be detectable 
as far away at low latitudes as it would be at closer to 
the Galactic pole. For example, an object with r^ « 30 



pc and My « -3.9 at 6 = 74° can be detected with 90% 
efficiency as far as ~ 120 kpc, while the same object at 
b — 31° has a 90% efficiency at ~ 95 kpc. 

To anticipate the effect that varying Galactic fore- 
ground will affect future dwarf searches in data that goes 
closer to the Galactic plane than SDSS, we also repeat 
the simulation of My ~ —3.9 galaxies at a foreground 
density of 10,000 stars per square degree, approximating 
a latitude of ~ 15°. The 50% and 90% detection effi- 
ciency contours (red) for these simulations are also shown 
on Figure [17] This further reduces the 90% detection 
distance of our example object to ~ 80 kpc but demon- 
strates that future surveys should still detect dwarfs at 
relatively low Galactic latitudes, barring extinction ef- 
fects. 

To check that latitude has the same lack of effect over 
different magnitudes, we calculate the 50%, 90%, and 
99% detection efficiency distance for a r/i = 100 pc object 
over the magnitude range —1.5 < My < —7.0 at b — 
31°, 53°, and 71° (Figure [T5] note that in this figure 
the distance scale is linear). The 53° and 71° curves 
are indistinguishable and the b = 31° curve is typically 
less than 20 kpc lower. This demonstrates that over the 
DR6 footprint, latitude does not play an important role 
on average in the detectability of objects. However, if we 
are unlucky, then individual objects could by chance lie 
in directions of unusually high foreground counts. 

6.5. Comparing Rgg, Rgo and R^q 

We use the analytical expression derived in §6.4 to esti- 
mate the distance at which each of the Milky Way dwarfs 
would be detected with 50% efficiency, R^o (Figure UM 
grey dots). Because R^q depends on half-light radius, as 
well as luminosity, we also for reference show R^q for ob- 
jects with rh = 250 pc (red) and rh = 50 pc (blue). These 
lines show that the detectability of the lowest luminosity 
dwarfs is severely reduced for large scale sizes. Objects 
with Segue 1, Bootes II, or Willman 1-like luminosities 
would not have been detectable with scale sizes of 100 pc 
or larger, even at very nearby distances. This size bias is 
important to bear in mind, particu larly given that the 
three M31 satellites discovered by iMcConnachie et al] 
(|2008[ ) that highlight regions of dwarf galaxy parame- 
ter space that have not previously been observationally 
seen. 

Assuming the size-luminosity distribution of known 
satellites is representative of all satellites we can ignore 
size and approximate R^q as a function of My by linear 
fit to the i?5o of actual dwarfs: 

logi?5o = -0.187A/y + 1.420. 



For comparison, the iKoposov et al.l (|2007aD equivalent 
of i?5o is shown in Figure [TO] (dotted line) obtained from 
Table 3 in their paper, and discussed further in S}7] For 
reference, the actual distances to the MW dwarfs are 
shown as black dots. 

A more useful quantity might be Rcompiete, the max- 
imum distance at which objects can be detected. Al- 
though we choose to use "complete" to refer to 90% de- 
tectability, complete could be defined as, say 90% or 99% 
efficiency. This can also be approximated by a linear fit 
to the results for the actual MW dwarfs: 



logi?9 



-0.204Afy + 1.164 
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Fig. 16. — Left: My = —3.9 panel of Figure [TSl binned to higher resolution. Center: Same as left panel but using the analytical 
expression to estimate detection efficiency. Right: Residual of the model-analytical efficiency. 




og(Distance) (kpc) 

Fig. 17. — The same as the My = -3.9 panel of Figure [TSl but 
showing the 50% and 90% contours of the simulations at 6 = 74° 
(blue), 31° (orange), and ~ 15° (red). 

or 

log-Rgg = -0.217Afy + 1.005. 

These relationships again assume that the known satel- 
lites are typical of the MW satellite population as a 
whole, i.e. objects with luminosities comparable to Segue 
1 are similar in size to Segue 1 and not significantly 
larger. 

Returning to Figure [TSl we see the distance range over 
which detectability changes from 99% to 50%. Rgo is 
typically ^ 20 kpc closer than R^q, and Rgg is ~ 20 kpc 
closer still. From this figure we can also see that objects 
brighter than My ~ —6.5 mag {My ~ —5.9 mag) are 
detected with 99% (90%) efficiency out to 300 kpc. We 
can infer from this that all dwarfs within the MW virial 




Fig. 18. — 50%, 90%, and 99% completeness distances for a 
rj, = 100 pc object as a function of magnitude at three Galactic 
latitudes: 31°, 53° and 71°. The 53° and 71° curves are virtu- 
ally indistinguishable showing that latitude does not significantly 
impact satellite detection over latitude ranges of DR6. 

radius brighter than My « —6.5 mag are known, and 
any satellites still undetected are likely to be reminiscent 
of objects like Coma Berenices, Bootes II or Segue 1 at 
distances greater than ~ 40 kpc. An ultra-faint satellite 
such as Segue 1 can only be detected with 50% efficiency 
out to ~ 40 kpc; there may be many more such objects 
beyond this distance. 

The code for interpolating detection efficiency from 
our simulations as well as the analytical function will be 
made available for download at the Astronomical Jour- 
nal website. Interpolation will give more accurate results, 
but the analytical function will provide flexibility for cus- 
tomization to suit individual needs and implementation 
of any future improvements. 

6.6. Caveats 
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Fig. 19. — Comparison of the 50% detection distance as a func- 
tion of magnitude for K07 (dotted) and for our analytical efficiency 
using r^ = 250 pc (red) and r/^ = 50 pc (blue). MW dwarfs are 
shown as filled circles. 

An underlying assumption of our simulations is that 
the DR6 point source catalog is uniformly 100% complete 
to r = 22.0. This assumption may result in optimistic 
detection efficiency estimates for the faintest and furthest 
systems. These faintest and more distant systems would 
also be subject to the human element; a real object may 
detected by the algorithm but on visual inspection be dis- 
regarded as background galaxy cluster or other contam- 
inant. Finally, sources in our sim ulated dSphs are dis - 
tributed circularly symmetrically. iMartin et al.l (120081 ) 
find that the ultra-faint satellites are in fact quite ellipti- 
cal, which due to our circular Plummer smoothing kernel, 
possibly results in overestimated efficiencies for objects 
like Hercules, UMa and UMa II with ellipticities of 0.68, 
0.80 and 0.63 respectively. 

7. COMPARISON WITH KOPOSOV ET AL. (2008) 

Besides IWillman et al.l ()2002l) from which this work fol- 
lows, two other surveys have recen tly uniformly sea rched 
SPSS for MW satellit es, namely iLiueFaLl (|2008[) and 
iKoposov et~an (|2007af ). iLiu et all (|2008[ ) conducted a 
straig htforward search and pr esented five satellite candi- 
dates. iKoposov et al.l ()2007a[ ) present a study compara- 
ble to this work, a nd follows from the "Field of Streams" 
(jBelokurov et al.l [2006a) that led to the discoveries of 
nine of the new Milky Way satellites. Here we compare 
our work in detail with Koposov ct al. (2007a) and sum- 
marize the m ain differences i n Tabl e [H 

The aim of IKoposov et al.l (|2007af ) was to present a lu- 
minosity function of the MW satellites corrected for lumi- 
nosity bias. Their analysis discovered two new extremely 
faint globular clusters, Koposov 1 and 2 (Koposo v et al.l 
2007b[). In principle, our analysis is quite similar to 



Koposov et al.l (|2007af l. henceforth K07, in that they ap- 
ply a color cut, smooth the stellar counts and look for 
statistically significant overdensities. There are several 
distinctions however that we detail below. 

K07 employed a g — r < 1.2 color cut to remove 
a substantial fraction of MW foreground stars and a 
r < 22.5 cut to limit the influence of background galaxies 



and increasing uncertainties/incompleteness. Our color- 
magnitude cuts are tailored to old stellar populations 
at 16 different distances which serve to eliminate more 
foreground stars than the looser K07 cut. The looser 
K07 color cut leaves enough stars that a complicated set 
of detection thresholds is unnecessary, whereas we must 
consider the effects of non-Gaussianity in low densities 
(see i j3.4[) . K07 deals with background galaxy clusters, a 
major source of contaminant overdensities, by producing 
a galaxy clustering significance in the same manner as 
the stellar clustering; anywhere that a stellar overden- 
sity occurs without a corresponding galaxy overdensity 
is much more likely to be a true stellar overdensity. Our 
algorithm only includes stars as faint as r = 22.0 and 
as such we have fewer galaxy cluster contaminant detec- 
tions. 

The most substantial difference between our work and 
K07 is how we derive the detection limits of our algo- 
rithms. Like our work, K07 simulated artificial galaxies 
to explore the detection efficiency as a function of size, 
distance and luminosity. K07 simulated 8,000 galaxies 
over a similar range of parameters as our study, but with 
only ~ 8 objects per 0.3 log{d) x 0.3 \og{rh) x 0.8 mag 
bin. There is considerable noise evident in the detection 
limits (their Figure 6), and all of the new satellites ap- 
pear to lie on the edge of detectability. K07 observed 
a steep, but finite, transition from unity to zero detec- 
tion efficiency which they attributed to the large range of 
distances that fall within each size-luminosity bin. How- 
ever, as discussed in our §5.21 a large number of simu- 
lated galaxies are essential for each permutation of dwarf 
galaxy parameters to effectively map their detectability. 
Our high resolution detection maps {^ 500 objects per 
0.1 log{d) X 0.1 log(r;i) X 0.8 mag bin, x3 latitudes) 
show that the detectability of a dwarf drops off slowly 
with size and distance, and that only Leo T, Leo IV, 
Leo V, Boo II and Willman 1 lie close to the edge of 
detectability. The difference between 90% and 10% ef- 
ficiency typically occurs over 0.2 dex in distance (kpc) 
and 0.3 dex in size (pc; see Figure [T5|l . 

Both the K07 and our detection limit calculations suf- 
fer from the implicit assumptions that the SDSS point 
source catalog is complete to the photometry limit and 
that dwarfs are circularly symmetric. These two as- 
sumptions yield detection limits that may be optimistic. 
The K07 study includes stars to a limiting magnitude of 
r — 22.5 mag, a half magnitude fainter than our limit of 
r = 22.0 mag. We thus expect that the completeness as- 
sumption may impact their calculated limits more than 
ours. 

To directly compare the effectiveness of both algo- 
rithms we return to Figure 1191 showing the distance at 
which an object is detected with 50% efficiency, R50 as a 
function of magnitude. The dashed line shows the K07 
i?50 which was determined by fitting a limiting magni- 
tude and surface brightness to the seven distance pan- 
els in their Figure 10. The dots show our R^q derived 
from the analytical efficiency function for each of the 
MW dwarfs in DR6, while the red and blue curves show 
i?50 calculated for objects of rn = 250 pc and r/j = 50 pc 
respectively. 

Although this comparison roughly shows that we have 
comparable limits, our calculated detection efficiencies of 
each dwarf are all greater than 90% while Table 2 of K07 
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lists efHciencies as low as 47% (neglecting Boo II) . While 
Bootes II is not detected with the standard algorithm of 
K07, it is a comparatively strong detection in our algo- 
rithm. We note that the tabulated K07 efficiencies for 
the known MW dwarfs appear inconsistent with their 
fitted i?5o, which places some dwarfs much closer than 
i?50 than their actual efficiencies would indicate. The 
increased dwarf detectability of our survey is owing to a 
combination of different techniques and our less stringent 
detection threshold. We set our thresholds to strictly 
eliminate truly random false positives expected while still 
yielding new candidates, and hence have ~ 30 unknown 
detections above our thresholds. Although upon visual 
inspection of their CMDs many of these detections ap- 
pear unlikely to be new dSph satellites, they may also be 
tidal debris or distant galaxy clusters. However, K07 set 
their detection thresholds just loose enough to retain all 
known objects; UMa is their weakest detection and there 
are only three unknown detections above this threshold. 

8. THE STILL MISSING SATELLITES 

A substantial driving force of this work is the Missing 
Satellite problem, which the discovery of so many new 
objects in the space of three years has shown is far from 
observationally exhausted. There are still large regions of 
parameter space where objects are undetectable, so there 
can easily exist more objects within the DR6 coverage 
that remain hidden. Future surveys such as the Stromlo 
Missing Satellites (SMS) Survey and Pan-STARRS may 
be able to detect some of these objects, and we can use 
our model of detectability to estimate how many there 
may be. 

We use a simplified version of the approach used in 
iTollerud et all ^QM). We first assume that the ra- 
dial distribution of dwarf galaxies matches that of all 
well-resolved subhalo es of the Via L a ctea simulation 
(jDiemandet al.ll2007D . iTollerud et all (|2008D discusses 
this assumption in detail; we realize this may not re- 
flect the true MW dwarf distribution, but our qualita- 
tive results are fairly robust to the assumed profile. For 
each satellite detected in DR6, we then determine Rgo 
(or -R99), the maximum distance to which a satellite of 
similar properties would be detected with 90% (99%) ef- 
ficiency. For each value of Rgo {R99), we determine from 
the Via Lactea subhalo radial profile what fraction of 
satellites should be within this distance, and weight each 
satellite accordingly. 

Using all dark matter subhaloes with more than 1000 
particles within r„ir = 289 kpc, arid ado pting a MW 
virial radius of 258 kpc (jKlvpin et al.ll2002i) . we estimate 
~ 13 (~ 24) satellites within the MW virial radius in 
the DR6 footprint. Twelve of these would be the known 
objects Bootes, Draco, Canes Venatici I and II, Coma 
Berenices, Leo I, Leo II, Leo IV, Leo V, Hercules, Ursa 
Major, and Ursa Major II, leaving one (12) possible miss- 
ing satellite(s). From our simulations, objects brighter 
than My — —6.5 mag are detectable with > 99% effi- 
ciency out to the virial radius, so we would expect that a 
relatively small number of the faintest systems are miss- 
ing. These missing satellites may be amongst our candi- 
dates, or be Coma Berenices like or fainter objects in the 
outer halo. Whether or not future searches reveal such 
objects may validate the assumed radial distribution. If 
we assume an isotropic sky distribution of satellites, ~ 13 



(~ 24) objects within DR6 equates to ^ 52 (~ 96) across 
the whole sky. 

If we include the ambiguous objects Segue 1, Willman 
1, and Bootes II 1 in the calculation, then the -R90 (Rgg) 
DR6 estimate would be ~ 56 (^ 85) satellites only 15 of 
which are known, or ^ 224 {^ 340) across the sky. The 
ambiguity of Segue 1, Willman 1, and Bootes II has con- 
siderable effect on the extrapolated MW census, under- 
scoring the need for an understanding of these extremely 
faint systems. 

These estimates assume that the sizes and luminosi- 
ties of the known satellites in DR6 are representative of 
the the MW satellite population as a whole. Based on 
our detection limits we cannot make any statements re- 
garding extremely diffuse, low luminosity systems that 
are undetectable by SDSS. The results also depend on 
the radial distribution assumed. If we instead assume 
that the Milky Way's dwarf population follows the ra- 
dial distribution of the MW dSphs known prior to 2004, 
then our Rgo inferred total number of dwarfs (with size- 
luminosities similar to known) within DR6 is 12, or 25 
including Segue 1, Willman 1 and Bootes II. This implies 
that all or most satellites within DR6 would be known. 

To Uerud et al.l (120 081 use the detection limits of 
Kopo sov et al.l (|2007al) to similarly estimate the true 
number of satellites within DR5 for a number of scenar- 
ios. The most comparable scenario to our assumptions 
(a limiting distance of 300 kpc, including all satellites 
except Segue 1 which is not in DR5) gives a result of 
322^75 satellites, consistent with our results of 232 for 
Rgo and 344 for Rgg. 

9. ALLOWING A LITTLE LATITUDE 

Substantial effort has gone into the observation and in- 
terpretation of the spatial distribution of the satellites of 
disk galaxies, in particular that of the Milky Way satel- 
lites. However, there is neither agreement on whether 
the Milky Way satellites have a truly anisotropic spa- 
tial di stribution, nor whether we expect them to. Pre- 
SDSS, iKroupa et all (|2005D found that the distribution 
of known MW satellites could be described by a disk of 
finite width, aligned almost perpendicularly to the MW 
disk. This was in a greement with the "Holmberg" ef- 
fect (|Holmberg|ll969f ) , that the closest satellites to a host 
galaxy were observed to be preferentially aligned with the 
minor axis of the host. This disk- like distribution seemed 
incompatible with ACDM, but Ka ng et al.l (|2005( ) rea- 
soned that if satellites follow the distribution of the host 
dark matter profile rather than that of the substructure 
then the dozen observed MW satellites could statisti- 
cally lie in a disk-like st ructure, al though the orientation 
of this disk is arbitrary. iPiatek et al.l (|2007f ) used proper 
motions derived from HST observations to show that this 
"Great Disk of Milky Way Satellites" was not a persis- 
tent structure; the orb its of the d\yarfs w ould not contain 
them within this disk. iMetz et al] (j2008f ) refute this con- 
clusion, finding instead that the orbital poles of most 
MW satellites place them in a rotationally supported 
disk-of-satellites. Studies on the satellites of other galax- 
ies fro m SDSS also yield confiicting results. 'Bailin et aU 
(|2007t l affirm the Holmberg effect while Brainerd (200^ 
find that satellites lie prefer entially along th e major, not 
minor, axis of the host. .Zentner et al.l (|2005f ) re-examine 
the problem from a theoretical point of view, stating that 
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TABLE 4 

Summary of Comparison with Koposov et al. (2008). 



This work 



Koposov et al. (2007) 



Survey Area 

Source cuts 

Smoothing kernel 

Threshold 

Modeled detection limits 

Number of Simulations 



DR6 - 9,500 sq deg 

Isochrone template at 16 dist intervals 

4.5' Plummer profile 

Multiple, function of foreground density 

HST obs of 3 MW dSphs 

3,825,000 



Efficiency map bin size (log(r(j)x 0.1x0.1x0.8 

log(d) X Mv) 

Simulation density (n per 0.3 X 0.3 X 0.8 ~ 4500 (x3 latitudes) 

bin) 



DR5 - 8,000 sq deg 

g-r <1.2 and r < 22.5 

o- = 2', 4', 8' - cr = 60' Gaussians 

Fixed, considers background galaxies 

M92 locus 

8,000 for general simulation + 1,000 each 

for known dwarfs within DR5 

0.3 x 0.3 X 0.8 



DM substructure is not completely isotropic and that 
the MW satellite distribution can, albeit with a very low 
probability, be drawn from a DM subhalo distribution. 

Within the standard ACDM structure formation sce- 
nario, satellite galaxies without dark matter could be 
formed in gas-rich tidal tails during vigorous earl y 
galaxy-galaxy interactions (jOkazaki &: Taniguchill200(i| ). 
Families of such tidal dwarfs would have correlated or- 
bital angular momenta and may appear as disk like ar- 
rangements about some hosts. This would support the 
apparent disk-of-satellites (Kroupa et al. 2005; Metz et 
al. 2007) and its correlated orbital angular momenta 
(Metz et al. 2008). It is therefore crucial to further con- 
strain the spatial and orbital angular momentum prop- 
erties of the satellites to reveal their true nature which 
is intimately related to the formation of the MW. 

A caveat of past studies of the MW dwarf distribution 
is that the sky had not been uniformly searched for satel- 
lites and the effect of Galactic latitude on the observabil- 
ity of dwarfs h ad not been thorough ly quantified. In their 
detailed study, iKlevna et al.l ()1997t ) showed that latitude 
strongly affected the detectability of Milky Way satel- 
lites with their technique. Our uniform study of SDSS 
DR6 takes their approach a step further and provides 
a detailed quantitative description of dwarf detectabil- 
ity over the footprint of our survey. We established in 
i i6.4l that the average detectability of the known satellites 
does not significantly vary over the DR6 footprint. We 
can thus compare the latitude and longitude distribution 
of the Milky Way satellites within the DR6 footprint with 
that expected if they are randomly distributed. We per- 
form a Kolmogorov-Smirnoff test to determine whether 
the satellites detected in DR6 show statistically signif- 
icant spatial anisotropy. Figure [20] shows the cumula- 
tive distribution (black lines) of latitude (top) and longi- 
tude (bottom) by area of the DR6 footprint. Overplot- 
ted on both panels are the cumulative distributions of 
MW satellites, both ignoring the ambiguous objects Boo 
II, Willman 1 and Segue 1 (blue) and including them 
(red). A KS test on these distributions with the entire 
DR6 area yields a probability of isotropic distribution of 
0.16 over longitude and 0.72 over latitude, or 0.06 and 
0.79 if we include Segue 1, Willman 1, and Bootes II. 
We also randomly pick 12 (or 15) points from DR6 cov- 
erage, weighted by area, and repeat the KS test 1,000 
times. The mean resulting probabilities of isotropic dis- 
tributions are 0.45 ± 0.29 over longitude and 0.62 ± 0.27 
over latitude, or 0.34 ±0.28 and 0.55 ±0.28 with Segue 1, 




(deg) 

Fig. 20. — Top: Cumulative histogram of the Galactic latitude 
of the DR6 footprint weighted by area (black). Blue shows the 
cumulative histogram of the latitudes of confirmed dSphs within 
DR6 and the MW virial radius (Boo, Dra, CVn, CVn II, ComBer, 
Leo I, Leo II, Leo IV, Her, UMa, and UMa II). Red shows the same 
but including Boo II, Willman 1, and Segue 1. Bottom: Same as 
top panel, but for Galactic longitude. 

Willman 1, and Bootes II. Hence, a conclusive result on 
the isotropy of the MW satellites awaits further data as 
our test shows that considering DR6 alone either scenario 
is plausible. 

10. TOWARDS SKYMAPPER 

SDSS Data Release 6 covers ^ tt steradians of Northern 
sky near the North Galactic Pole, and at least nine new 
dSph companions have been found in this area. The new 
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ANU SkyMapper telescope will survey the entire south- 
ern sky to a similar depth as SDSS over the next five 
years, so we can naively expect to find around twenty- 
five new Milky Way satellites. With a detailed, system- 
atic search covering around three quarters of the sky, we 
will for the first time be able to conclusively compare 
the MW satellite galaxy population with theoretical pre- 
dictions. The apparent anisotropy of the satellites will 
be conclusively confirmed or ruled out and we will con- 
tinue to discover the most dark matter dominated stellar 
systems nature produced. 

Designed as a replacement for the Great Melbourne 
Telescope that was destroyed by the Canberra bush fires 
in January 2003, SkyMapper is a new wide-field auto- 
mated survey telescope currently being built at Siding 
Spring Observatory, NSW. It will feature a 1.33m pri- 
mary mirror and a 0.69ni secondary with an effective 
aperture of 1.13m and a 5.7 square degree field of view. 
The imager consists of a 4 x 8 mosaic of 2048 x 4096 
pixel CCDs (jKeller et al.l 120071 ). SkyMapper's primary 
purpose is the Southern Sky Survey which will be used 
to study objects from Trans-Neptunian Objects to high 
redshift quasars. It's five year mission: to survey the en- 
ti re 27r steradi a ns of t he southern sky at multiple epochs. 

iKeller et al.l ()2007f ) states an average seeing at Sid- 
ing Spring of ~ 1.5 arcsec as derived from observing 
logs of the Anglo- Australian Telescope (AAT), on par 
with the median 1.4 arcsec se eing of the SDSS site 
(jAdelman-McCarthv et~aLl 120071) . The sky will be ob- 
served at six epochs, and at completion 90% of the 
sky will be observed at least five times, and 100% ob- 
served at least three times. The expected survey depth 
(S/N — 5, t=110s per epoch) is given in Table [H 
Also include d are the corresponding magni tude limits 
from SDSS (|Adelman-McCarthv et all 12001 . The sur- 
vey aims to provide astrometry to better than 50 mas, 
as compared with 100 mas for SDSS. Further information 
is available at the SkyMapper website^. 

11. CONCLUSION 

The dwarf galaxy satellites of the Milky Way provide 
excellent opportunities to further our understanding of 
galaxy formation and near-field cosmology. They can be 
resolved into individual stars allowing detailed studies of 
their structure, kinematics and composition. They have 
also been cause for concern regarding their interpretation 
as the luminous components of dark matter substructure; 
it has been argued that the number and spatial distri- 

^ http://www.mso.anu.edu.au/SkyMapper 



bution of these satellites are inconsistent with ACDM 
structure formation scenarios. The commencement of 
the Sloan Digital Sky Survey triggered a cascade of dis- 
coveries, with fourteen new satellites discovered. The 
limited spatial coverage and photometric depth of SDSS 
suggests that many, if not most, MW satellites are still 
yet to be discovered. The coming years are likely to 
bring the MW satellite census towards completeness as 
new survey telescopes such as SkyMapper, Pan-STARRS 
and LSST come online. 

We present here the method used to search the Sloan 
Digital Sky Survey Data Release 6 for ultra-faint Milky 
Way satellite galaxies. By screening for stars consistent 
with an old population at a fixed distance, we enhance 
the signal of a dSph over the Milky Way foreground. 
Smoothing with a kernel corresponding to the expected 
surface density profile further elevates the dSph above 
the foreground, and our comprehensive thresholds ac- 
count for varying stellar density and more diffuse objects. 

Applying our algorithm to SDSS DR6, we recover the 
"classical" and recently discovered dSphs, as well as 17 
globular clusters and two open clusters. We also have 
30 unidentified detections, some of which may be new 
satellites. The discovery of Leo V demonstrates the dif- 
ficulty in following-up dwarf candidates; while we detect 
Leo V, there are several unknown detections of greater 
significance that may prove to be something. However 
observing these weakest can didates is a rather hit and 
miss affair, as pointed out bv lBelokurov et all (j2008[ ). 

To compare the known dwarf galaxy population of 
the Milky Way with predictions, its essential to have a 
very well-defined dwarf selection function. To do this, 
we thoroughly model the detection efficiency of systems 
covering a wide range in parameter space by simulating 
more than 3,000,000 galaxies. We fit various functions 
to the resulting detection efficiency contours to semi- 
analytically describe efficiency as a function of magni- 
tude, size, distance and Galactic latitude. Using the re- 
sults of our detailed investigation of dwarf detectability, 
we show that: 

• Assuming a Via Lactea subhalo radial distribution 
and that Rcompiete — R90, there should be '^ 13 
satellites with DR6, 12 of which are known. If we 
include Segue 1, Willman 1, and Bootes II in this 
calculation, this estimate jumps to ~ 56, only 15 
of which are known. 

• Dwarf detectability shows a smooth transition from 
100% to 0% over size and distance. For example, 
the distance at which a CVn Il-like object is de- 
tected with 90% efficiency is 200 kpc, compared to 
316 kpc for 10% efficiency. 

• Galactic latitude does not significantly impact the 
detection of satellites over the DR6 footprint, and 
surveys of similar quality should still detect dwarfs 
as low as 6 « 15°. All of the satellites discovered 
in SDSS would have been detected at any latitude. 

• The census of MW satellites brighter than My — 
—6.5 should be complete out to 300 kpc, and all 
objects brighter than My = —5 would be detected 
with at least 50% efficiency out to this distance. 
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• Given the present data, the spatial anisotropy of 
the MW satellites within DR6 cannot be confirmed 
or ruled out. 

We provide several different parameterizations of our 
detection limits to facilitate comparisons between the 
known Milky Way dwarf galaxy population and predic- 
tions. We provide software that returns the detection ef- 
ficiency of a dwarf galaxy as a function of its luminosity, 
scale size, distance, and latitude. There are two different 
codes provided for this; one is based on an analytic de- 
scription of our detection limits and the other provides a 
direct interpolation from our 3,825,000 simulated galax- 
ies. We also provide a linear fit as a function of My of 
the distance out to which dwarfs are detected with each 
of 50%, 90%, and 99% efficiency. These fits assume an 
underlying dwarf galaxy population with combinations 
of sizes and luminosities similar to those known. 

2009 will bring about the beginning of the South- 
ern Sky Survey, and with it a way to uniformly search 
a further 20,000 square degrees of sky for new MW 
dwarfs. We will apply our algorithm to the incoming 
data and produce the most complete and well charac- 
terized census of the MW neighbourhood possible to 



date. We may also implement improvements to the al- 
gorithm to o ptimize for stellar streams or young stellar 
populations. iMcConnachie et all (|2008f ) have shown that 
M31 satellites occupy as yet unexplored size- luminosity 
space around the Milky Way. Surveys beyond SDSS and 
SkyMapper, such as Pan-STARRS and LSST, wiU be 
needed to carefully search for such systems. Even with 
our carefully characterized detection limits the true num- 
ber of MW satellites remains highly uncertain. 
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